Overdamped dynamics of long-range systems on a one-dimensional lattice: 
Dominance of the mean-field mode and phase transition 



c/3 



Shamik Gupta 1 Alessandro Campa 2 [j] and Stefano Ruffo 1 '^ 
1 Laboratoire de Physique de I'Ecole Normale Superieure de Lyon, 
Universite de Lyon, CNRS, 46 Allee d'ltalie, 69364 Lyon cedex 07, France 
2 Complex Systems and Theoretical Physics Unit, 
Health and Technology Department, Istituto Superiore di Sanita, and INFN Romal, 
Gruppo Collegato Sanita, Viale Regina Elena 299, 00161 Roma, Italy 
3 Dipartimento di Energetica "Sergio Stecco" and CSDC, Universitd di Firenze, 
CNISM and INFN, via S. Marta 3, 50139 Firenze, Italy 
(Dated: October 1, 2012) 

We consider the overdamped dynamics of a paradigmatic long-range system of particles residing 
on the sites of a one-dimensional lattice, in the presence of thermal noise. The internal degree of 
freedom of each particle is a periodic variable which is coupled to those of other particles with an 
attractive Xy-like interaction. The coupling strength decays with the interparticle separation r 
in space as l/r Q ; < a < 1. We study the dynamics of the model in the continuum limit by 
considering the Fokker-Planck equation for the evolution of the spatial density of particles. We 
show that the model has a linearly stable stationary state which is always uniform in space, being 
non-uniform in the internal degrees below a critical temperature T = 1/2 and uniform above, with 
a phase transition between the two states at T = 1/2. Thus, the stationary state of the model is the 
same as that of its mean-field counterpart which has a = 0. We justify this mean-field dominance 
by performing linear stability analysis of both the uniform and non-uniform stationary solutions of 
the Fokker-Planck equation. Our analysis also allows us to compute the growth and decay rates 
of spatial Fourier modes of density fluctuations. These rates compare very well with numerical 
simulations. 
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I. INTRODUCTION 

In long-range interacting systems, the interparticle po- 
tential in d dimensions decays at large separation, r, as 
l/r", with a < d \u^- Examples are gravitational sys- 
tems [H , plasmas two-dimensional hydrodynamics Q , 
charged and dipolar systems Q, etc. Long-range sys- 
tems arc non-additive, and have equilibrium properties 
unusual for short-range systems, e.g., inequivalence of 
statistical ensembles 0. Besides, they show intriguing 
dynamical features, e.g., broken ergodicity [{J and slow 
relaxation towards equilibrium 0, [T3| ■ 

A paradigmatic example of long-range systems is the 
so-called Hamiltonian mean-field (HMF) model compris- 
ing a system of N interacting particles evolving un- 
der Hamilton dynamics (Tlj . The ith particle, i = 
1, 2, . . . , N, has an internal degree of freedom 9[, a pe- 
riodic variable of period 2ir. The Oi's are coupled to 
one another via an attractive Ay-like interaction ~ 
— cos(9i — 0j) between the ith and jth particles, with 
the coupling strength being of the "mean-field" type: 
It is equal for every pair of particles. As a result, the 
model is defined without requiring any lattice structure 
for the particles to reside on. In this model, a wide class 
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of initial distributions relaxes to equilibrium over times 
that diverge with the system size [10[ • As a consequence, 
the system in the thermodynamic limit does not ever at- 
tain the Boltzmann-Gibbs equilibrium, instead it remains 
trapped in intermediate quasistationary states. 

To explore how far the results obtained within the 
HMF model extend beyond its setting of mean-field 
coupling, the model has been considered on a one- 
dimensional periodic lattice on the sites of which the par- 
ticles reside. Particles on different sites are coupled as in 
the HMF model, with the difference that the coupling 
strength decays with their separation r along the lattice 
length as l/r a , with < a < 1. The resulting model is 
known as the cv-HMF model [l2[ . Within a canonical en- 
semble, this model has the same equilibrium properties 
as those of the HMF model [Hi]. This raises an immedi- 
ate question: Does this equivalence extend to dynamical 
properties as well? This work is a step towards answering 
this question. 

To address the issue, let us consider the simple setting 
of the cv-HMF model in the presence of thermal noise and 
dissipation. The ith particle has the equations of motion 



dBj_ 
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where pi denotes the momentum of the ith particle, while 
7 is the damping constant. The second term on the right 
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hand side of Eq. (|2|) is the force on the ith particle arising 
from the interaction potential: The quantity \j— i\ c is the 
closest distance on the periodic lattice between the ith 
and jth sites: 
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while 



i\ r = mm 



JV 



(\j-ilN-\j-i\), 



Vi. 



(3) 



(4) 



From Eq. ([5]), since a < 1, the cumulative interaction of 
one particle with all the remaining particles with aligned 
#'s would diverge in the limit N —> oo in the absence of 
the normalization N, which explains its inclusion [13j . In 
Eq. @, Tji(t) is a Gaussian white noise: 



(Vi(t)) = 0, 

(Vi(t)Vj(t'))=2T6 ij S(t-t') 



(5) 
(6) 



where T denotes the temperature, while the angular 
brackets denote averaging over noise realizations. In this 
work, we take the Boltzmann constant to be unity. 

Here, we will study the ovcrdampcd limit of the equa- 
tions of motion, (QJ and (J2J ; after a time rescaling, these 
equations then reduce to the following Langevin equation 
for the ith particle: 



N 



dtii _ 1 x - 
~dt ~ ~N ^ 
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+ Vi{t); a>0. (7) 



We analyze the dynamics (J7J in the continuum limit 
N —> oo, when the lattice becomes a continuous segment 
characterized by the spatial coordinate s e [0, 1]. In this 
limit, the local density of oscillators p(9; s, t) evolves in 
time following a Fokkcr-Planck equation. This equation 
allows a stationary state which is uniform both in 9 and 
s. By performing a linear stability analysis of such a 
state, we show that when it is unstable, different spatial 
Fourier modes of fluctuations have different critical tem- 
peratures below which the modes grow exponentially in 
time with different rates. The largest critical tempera- 
ture, T C fl = 1/2, corresponds to the zero Fourier mode 
or the mean-field mode which has no spatial dependence. 
Above this temperature, all the Fourier modes decay in 
time, thereby stabilizing the uniform state. Below T Cj o, 
our numerical simulations starting from a uniform state 
show that, after an initial exponential growth, all the 
non-zero Fourier modes decay at long times to zero. By 
contrast, the mean- field mode grows and reaches a non- 
zero value, corresponding to a non-uniform stationary 
state, i.e. a state which is non-uniform in 9, but uniform 
in s. Such a mean-field dominance at long times has been 
observed in the a-HMF model evolving within a micro- 
canonical ensemble, i.e. equations ([T]) and ([2]) with 7 = 
[l4j | . In that work, this was just a numerical observation 
with no theoretical justification. Here, within the canoni- 
cal and overdamped dynamics, we justify this mean-field 



dominance by performing a linear stability analysis of 
the non- uniform state. We conclude that for any a in 
the range < a < 1, the linearly stable stationary state 
of the model ([7]) in the continuum limit is always uniform 
in space, being non-uniform in 8 below T Ci o and uniform 
above, with a phase transition between the two at T c .o- 

The paper is structured as follows. In section UH we 
introduce the continuum limit of the dynamics ([7J , and 
perform the linear stability analysis of the uniform sta- 
tionary state. In section Hm we present numerical simu- 
lations in support of the analysis, in particular, to show 
the agreement of the growth rates of the unstable modes, 
the decay at long times of the non-zero modes at all tem- 
peratures T < 1/2, and the associated dominance of the 
mean-field mode. In section IIV1 we perform the linear 
stability analysis of the non-uniform stationary state. We 
end the paper with conclusions and perspectives. 



II. CONTINUUM LIMIT: UNIFORM 
STATIONARY STATE AND DYNAMICS OF 
FOURIER MODES 

In the continuum limit N — > 00, the lattice of the 
system ([7]) is densely filled with sites. Let us introduce 
the variable s = i/N to denote the spatial coordinate 
along the lattice length, such that in the continuum limit 
it becomes the continuous variable s £ [0,1]. In this 
limit, we define a local density of particles p(8; s, t) such 
that of all particles located between s and s + ds at time 
t, the fraction p(9; s, t)d9ds have their degrees of freedom 
between 9 and 9 + d9. This density is non- negative, 2tt- 
pcriodic in 9, and obeys the normalization 



d9 p(9;s,t) = 1 V s. 



(8) 



In the continuum limit, the equation of motion is given 
by 



— ±L± = K (a) J dO'ds' > p{e'; s', t) + v (s, t), 



where 



re(a) 1 = 



s+1 ' 2 ds' 



s-l 



/2 (w-sur 1 



and 



\s — s\ c = min(|s — s|, 1 — \s — s\). 
Also, we have 

( v (s,t)} = 0, 

(r)(s, t)r)(s', t')) = 2TS(s - s')S(t - t'). 



(9) 



(10) 



(11) 



(12) 
(13) 
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The time evolution of p(9;s,t) follows the Fokkcr- 
Planck equation 



dp 

at 



+T 



d9 2 ' 



(14) 



Its stationary solution, obtained by setting the left hand 
side to zero, is 



Po(9; s) ex exp 



k{q) 



T 



de'ds 



, cos(W - 0) 



Po(0'; s' 



(15) 

Consider the particular stationary state which is uni- 
form both in 9 and in the spatial coordinate s, 



1 

Z7T 



(16) 



This state, although stationary, might be destabilized by 
the thermal fluctuations inherent in the dynamics. Let us 
then analyze the stability, in particular, linear stability 
of the state with respect to small thermal fluctuations 
Sp(9; s,t). To this end, we write 

p(e;s,t) = ^+6p(6;s,t); 5p(6;s,t) < 1. (17) 

Substituting Eq. (fT7|) into Eq. (fT4|) and keeping terms 
to linear order in dp, we find that bp satisfies 



dSp n(a) 
~~dt ~ 2tt 



, , cos(6>' - 9) . . rT1 d 2 5p 

d9'ds' ^-5 P {6';s',t) + T^- 
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Expressing Sp in terms of its Fourier series with respect 
to the periodic variable 9 as 



(19) 



k= — oo 



we find from Eq. ([18)) that only the modes k = ±1 are 
affected by the coupling between the particles, and that 
6p±i satisfies 



dSp ±1 
dt 



-T5p ±1 



K{a) f j , Sp ±1 (s',t) 

(W-S\c) a 



ds 



(20) 



One thus gets a set of equations for each position s, all 
coupled together by the second term on the right hand 
side of Eq. (gDJ). For k ^ ±1, one has 



9Sp ±k 
dt 



= -TSp ±k ; k^l. 



(21) 



Thus, 5p± k (s,t) for k ^ ±1 decays exponentially in time 
as exp(— Tt), so that these modes cannot destabilize the 
uniform state (JTGj) . 



Since we have a periodic lattice, to solve Eq. (|20|) . we 
next consider the spatial Fourier series of <5/9_(_ 1 (s,t): 



Sp ±1 (s,t)= £ S P±hm (t)e^ 



(22) 



On substituting in Eq. (|20[) . we get 



dt 



T K (a)A m (a) I - 



(23) 



where 



A m (a) 



s+l/2 i2nm{s'-s) 

d s' 7T-, i-r-- ( 24 ) 

.-1/2 (Is'-Sic)" 



Note that A m (a) = A_ m (a). It is known that 
A m (o;) > 0, and that it is a monotonically decreasing 
function of m when m is an integer; moreover, A m (a) — > 
0asm-> ±oo [3]. 

Equation (|23|) has the solution 



Sp±i. m (t) = cxp ( 



'K(a)K m (a) 



T)t 



6p ±1<m (0). (25) 



It then follows from our linear stability analysis that the 
mth spatial mode of fluctuations Sp ±l m (t) either (i) de- 
cays in time, which happens at temperatures such that 

rp K(q)A m (q) 

T < 

haviors is achieved at the critical temperature T c , m for 
the mth mode, given by 



, or (ii) grows in time at temperatures 
«(a)A m (a) ^ rp ne |j 0r( j er i me between these two be- 



K(a)A„ 



(26) 



The mth mode at T < T c m grows in time as exp[(T Cim — 
T)t\. With A m (a) = A_ m (a), the Fourier modes m and 
— m have the same critical temperature. Since Ai m i(a) 
decreases on increasing |m|, we conclude that 



T c ,o > T c i > T C 2 > ■ ■ 



(27) 



Note that T Cj o is the critical temperature for the m = 
mode, i.e., the "mean- field" mode; we have 



«;(a)Ao(a) 1 
2 _ 2' 



(28) 



Thus, above T = 1/2, when all the spatial modes decay in 
time, the uniform state (fT6)) is stable. On the other hand, 
the state is unstable when T < 1/2. At temperatures 
between 1/2 and T Cj i, the m = mode grows in time, 
while all the other modes decay in time. At temperatures 
between T Cj 2 and T Cj i, the m = 0, 1 modes grow in time, 
while all the other modes decay in time, and so on. In 
Table HI we show representative numbers for these critical 
temperatures for a = 0.25, 0.5, 0.75. 
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a 


T c ,i 


T c , 2 




To, 4 


Tc,5 


0.25 


0.082555 


0.042070 


0.033719 


0.025764 


0.022664 


0.50 


0.186991 


0.122063 


0.103418 


0.087614 


0.079556 


0.75 


0.321783 


0.262301 


0.239974 


0.221807 


0.210691 



TABLE I: Critical temperatures T Cjm ; m = 1, . . . , 5 for a — 
0.25,0.50,0.75. Note that T c , = 1/2, independent of a. 
When a = 0, one has T c , m = for m =^ 0. When a ap- 
proaches unity from below, all the T c , m 's for m/0 approach 
T c , o ■ 



III. COMPARISON WITH NUMERICAL 
SIMULATIONS 

In order to test the theoretical predictions in the limit 
./V — > oo obtained in the preceding section regarding val- 
ues of critical temperatures, and growth and decay of 
Fourier modes depending on the temperature, we per- 
formed extensive numerical simulations of the dynamics 
for large but finite N. A standard procedure for simula- 
tions is to integrate the equation of motion (J7J for each 
of the N particles; this involves computing at every inte- 
gration step a sum over N terms for each particle, thus 
requiring a total computation time scaling with N as N 2 . 
In order to perform faster simulations, we adopted the al- 
gorithm discussed in Ref. fl5| to compute the sum in Eq. 
([?])■ This algorithm requires a computation time that 
scales with N as N\nN. The equations of motion were 
integrated using a fourth-order Runge-Kutta algorithm 
with time step equal to 0.01. The initial state of the sys- 
tem was chosen to be the uniform state (|16p , prepared by 
having each particle degree of freedom 9 uniformly dis- 
tributed in [— 7r, 7r], independently of the others. Thus, 
initially, the system is in a state which is uniform both 
in 9 and in space. We report here simulation results for 
N = 2 14 . 

In simulations, we monitor the observable 



i=i 



;m= 0,1,2,.. 



(29) 



in particular, ro(t) does not contain any spatial depen- 
dence, and, thus, characterizes the mean-field mode. In 
the continuum limit, we have 



.(*) 



dOds e^ e+2mns) p(9;s,t) . 



(30) 



Starting at time t = from the uniform state such 
that r m (0) = V m, and at a temperature T = 0.01 
with a = 0.25, we expect from Tabic Q] that the Fourier 
modes 0, 1, 2, 3 in particular are linearly unstable. Con- 
sequently, 7"o(i), ri(t), T2(t), J~3(i) are expected to grow in 
time. Indeed, Fig. QJa) displaying the simulation re- 
sults for r (t) , T\ (t) , r 2 (t) , r$ (t) does show that they all 
grow in time. However, in the long-time limit, we see 
that ro(t) saturates to a value very close to unity, while 



r i{t)i r 2{t),rj,(t) all decay to a value very close to zero. 
By repeating simulations at different temperatures T < 
1/2 for different values of a (see Fig. [TJb),(c)), we have 
confirmed that ro(t) always saturates in the long-time 
limit to a non-zero stationary state value rg* = rg*(T) 
that depends only on the temperature, while r m (£)'s, 
with m ^ 0, decay in the long-time limit to a station- 
ary state value very close to zero. These observations 
suggest that for temperatures T < 1/2 for all values of 
a, the stationary state is characterized by a distribution 
po which is non- uniform in 9, but uniform in s, and may 
be obtained from Eq. (fT5")) as: 



p (9) oc exp 
= exp 



f d9> ds > C -^ ] 
T J (|s -s|c) 

— (m x cos 9 + m y sin ( 



>(#') 



(31) 



where 



d9{cos9,sin9)po{9). (32) 



IV. NON-UNIFORM STATIONARY STATE AND 
DYNAMICS OF FOURIER MODES 

Our observations regarding saturation of ro(i) and de- 
cay of r m (£)'s, with m ^ 0, detailed in the preceding 
section, raise an immediate question: Can we show by 
performing a linear stability analysis of the state (|3Tj) at 
temperatures T < 1/2 that all the modes of fluctuations 
decay in the long-time limit to zero, so that the state is 
linearly stable for all values of a? 

Before proceeding with the stability analysis, we note 
that our system of study has full rotational invariance in 
9 at each space point, so that in the following, we take 
rriy = and < m x < 1 without any loss of generality; 
then, we have 



p {9) = A exp 
where A is the normalization: 
.4 



m x cos 9 



1 



(33) 



(34) 



Note that the quantity m x is determined sclf-consistcntly 
by combining Eq. (|3"3")l with Eq. ([52"]) : one gets 



h{m x /T) 



(35) 



I {m x /Ty 

where I n (x) is the modified Bessel function of order n: 



In(x) 



1 

2^ Jo 



d9 e xcos 9 cos n9. (36) 



We thus have A = [2tt I a {m x /T)] \ 
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To study the linear stability of the state (|33[) , we write 
p(e;s,t)=p (e) + 5p(9;s,t); 6p(9; s,t) « 1, (37) 
where we expand 6p(9; s, t) as 



5p{6;s,t)= J2 8 Pn (6,t)i 



il-wms 



(38) 



Using Eq. (|38]) in Eq. (fj"4|) . we find that Sp„.(6. t) satisfies 

dSp m d(sm6 8p m ) d 2 6p m 
= m x — h 1 



9 .( 



at 

where we have defined 

A m (a) = re(a)A m (a) 



69 ' ~ 69 2 
d9'sm{9' -6)5 Pm {9',t) p ),(39) 



(40) 



The eigenfrequencies associated with the linear equa- 
tion (|39]) can be studied by looking for solutions of the 
form 



8p m {e,t) = 5p m (6,uj)e i 



(41) 



Before studying the spectrum of the eigenfrequency w, 
it is instructive to analyze the explicit solution dp m (9,0) 
corresponding to the neutral mode w = 0. We find from 
Eq. J39]) that 5p m (9,0) satisfies 



m x sin 9 5p m — X m (a) (m y cos 9 — m x sin 9 ) po 
,86 p. 



+T- 
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where C is a constant independent of 9, and 
(fh x ,fh y ) — J d9(cos 9, sin 9)Sp m (9, 0). 

Equation (|42|) has the solution 

Sp m {0,O) = «Tp m (0,0) e -^( 1 - cose ) 
Cm r e 

T Jo 

AX m (a)fh y ^l cos8 . „ 
H ^ — -e smf 



(42) 



(43) 



T 

AX m (a)fh x s 



r 



e T 



cos 



(1 - cos 6>). (44) 



The periodicity condition 8p m {9 + 2n, 0) = Sp m (9, 0) im- 
plies that C = 0. Normalization of p(9; x, t) for all times 
implies, using Eq. (|3T|) and the fact that po(^) is nor- 
malized, that J Q * d9 Sp m (9, 0) = V m. Then, from Eq. 
(|44j) with C = 0, we get 



2jt 



Jo 



AX m (a)m x 



d8t 



\\-cos0). (45) 



Using the above equation and Eq. (|35|) in Eq. (|44)) . we 
arrive at the following expression for <5p m (0,O): 



AX m {a) my., 

T 6 



[fdj (cos — m x ) + m B sin 0] . 



(46) 

The quantities in Eq. (|46|) arc determined self- 

consistently. They cannot both be equal to zero, as would 
have happened if 5p m {9, 0) were a constant, but which is 
not possible as the integral of 6p m (9,0) from to 2tt is 
required to vanish. Using Eqs. (031), (gBJ), (O and ([55]), 
we get the following two equations: 



- ~ -WoQ h(m x /T) 
m x I Q (m x /T)' 



~ A m (a) 

= 7T7 



T 



'l-T- 



m~ 



(47) 
(48) 



where < T < 1/2, and m x satisfies Eq. (|33|. Equation 
(|4T| is obviously satisfied for m 9 = 0. If m y ^ 0, then, 
using Eq. (J35j, we see that Eq. (]47)) is an identity for 
m = 0, while it has no solution for m ^ 0, which follows 
from the property that X m (a) < 1 for m =/= 0. Let us now 
analyze Eq. (|4"8"1) . It is satisfied for rh x = 0. If fh x ^ 0, 
it gives 



v/A„» (1 - T) - T, 



(49) 



giving a relation between m x and T, which has to be sat- 
isfied together with Eq. (|35[) . To check whether the two 
equations are consistent, first note that X m {a) > is a 
monotonically decreasing function of m, with Xo(a) = 1 
and Ymim^aa X m (a) = 0. Hence, to check the consis- 



tency, we have to consider in Eq. (|49|) any value for 
X m (a) between and 1. Actually, when A m (a) < 1, 
Eq. (|49"1) gives a real value for m x only for T in the 
range < T < , which is smaller than the range 

< T < 1/2. In Fig. H we plot m x (T) as given by 
Eq. (]4^]l for X m (a) = 1, and as given by solving the im- 
plicit equation (J35j. We see that the two curves do not 
intersect. This also implies that there cannot be any in- 
tersection for X m (a) < 1, since the right hand side of Eq. 
(|49]) decreases on decreasing X m (a). Thus, we conclude 
that Eq. (J15J for < T < 1/2 is not solved by m x that 
satisfies Eq. (|35|1 . In summary, the two equations, (|47j) 
and l|48p. are not satisfied by m x given by Eq. (|3"S"]) . when 
m > 0. When m = 0, there is only one solution, namely, 
m x = 0, and m y any non-zero value. For this solution, 
on using Eqs. (|33[) and pt?]) . we have 



<Vo(M) = ^po(0)m„sin< 



(50) 



which is just a perturbation of po(9) that corresponds to 
a rotation in space of all the rotators of the system. 
It is therefore natural that it corresponds to a neutral 
mode. 
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On the basis of the above discussion, we conclude that 
there is no temperature T < 1/2 such that the mth mode 
of fluctuations, for any m, has a vanishing cigcnfrcqucncy, 
excepting for the trivial one for m = that corresponds 
to a state obtained by a rotation in 9 space of all the ro- 
tators of the system. Thus, the mth mode either grows 
or decays in time at all temperatures T < 1/2, corre- 
sponding to cigcnfrequencies oj which have respectively a 
negative or a positive imaginary part. 

In order to examine the entire spectrum of the eigenfre- 
quency u>, we use Eqs. (f3U)) and (|4ip . and the substitution 
fx = iu to get the following equation for 8p m (0,u>). 

~ dismJUp^) d 2 Sp m 
pdp m = m x — 1- T- 
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d ( 



-A m (a)— y / d6'sm(e'-e)5p m (6',uj) ft,). (51) 



Let us define 



(m("',mi n) )= / de(cosne,smne)6p m (e,oj), 



(52) 



4"'= I 'decode Po {6) = Yr xl ?TV 

J Io(m x /T) 



where we note that mjiP = m x . On multiplying Eq. ()51j) 
in turn by cosnO and sinn#, integrating over 9 from 
to 27r, noting that po in even in 9, that 5p m (2Tr,Lu) = 
<5p TO (0,a->), and that pq(2t:) = po(0), we arrive at the 
following system of equations for n > 1: 



pm^ 



+ l -\ m {a)nm^ (m^ - m<? + V), (53) 

^nm x (rh^-V - m y n+1] ^j - Tn 2 m y n ^ 

+ lx m (a)n^( mx ^ +m ^ (54) 



where fh x ^ = m y J> = 0, and m x '' = 1. The operators on 
the right-hand sides of the above equations being non- 
hermitian have in general both real and complex eigen- 
values. 

Combined with Eq. ([33]) , the eigenvalues of the system 
of equations, ([53]) and (f54|) . can be obtained numerically 
after truncating the equations to a finite value of n, say, 
n max . As argued before on the basis of properties of 
A m (a), in order to probe the behavior of a m ^ mode 
of fluctuations, it suffices to consider any real positive 
value < 1 for A m (a), while for the m — mode, we have 
A m (a) = 1. The results for the eigenvalue spectrum of 
Eqs. ([52]) and (JM]) are shown in Figs. [3] and Fig. 0] for 
a temperature T < 0.2 and at the temperature T = 1/2, 
respectively. 

Let us discuss the results displayed in Fig. |3J From 
the m/0 mode results, displayed in panels (a) and (b), 



we see that the eigenvalues of the x- and the y-system 
are both real and complex, but the important thing to 
note is that the real eigenvalues are all negative while the 
complex ones have strictly negative real parts. In order 
to demonstrate the former fact, a zoom into the region 
near the zero of the Re(p) axis, shown in the insets, illus- 
trates that the eigenvalue closest to of both the x- and 
the y-system has a negative real part, so that the remain- 
ing eigenvalues having larger real parts in magnitude are 
thus all negative. The insets also show that the eigen- 
value with the smallest negative real part converges in 
magnitude on increasing the truncation order n max . For 
higher n max , only the number of eigenvalues increases; 
in particular, only the number of real eigenvalues with 
larger magnitude increases, while there is still only one 
complex eigenvalue which has the largest negative real 
part. Computing the eigenvalues for other A m 's and tem- 
peratures T < 1/2, we find that the number of complex 
eigenvalues is always small (with the number increasing 
for smaller A m 's and T"s), and these eigenvalues always 
have large negative real parts. Since the mth mode of 
fluctuations has the time dependence ~ e M * (recall Eq. 
(|4"Tj) and the definition of p, = ioj), and since all values 
for p have strictly negative real parts, it follows that at 
temperatures T < 1/2, the mth mode of fluctuations, 
with m ^ 0, decays in the long-time limit to zero. 

Figure[3J panel (c) shows that the behavior for the m = 
mode is very similar to that observed in panels (a) and 
(b) for m / 0, discussed above. Panel (d) also has the 
same general behavior, excepting that now the eigenvalue 
closest to the origin of the Re(p) axis, shown in the inset 
as a function of n max , is very slightly (~ 10~ 7 ) positive. 
Actually, this eigenvalue should be exactly zero, being 
the one corresponding to rotational invariance that we 
discussed earlier (see the discussion following Eq. (|49p). 

In fact, using Eq. ([50)1 in Eq. (|52")) to compute 



■An) 

ly , 



and 

using it in Eq. (|54|) , it may be checked that its right hand 
side is zero for any n, that is, the particular solution (|50[) 
gives p = 0. The slight positive value is just a numerical 
artifact. 

We now discuss the results displayed in Fig. @] From 
the m/0 mode results, displayed in panels (a) and (b), 
we see that the eigenvalues of the x- and the y-system 
are only real and negative, with no complex eigenvalue. 
The eigenvalue closest to the origin of the Re(p) axis, 
shown in the insets as a function of n max , dominates the 
behavior of the mth mode of fluctuations in time, thereby 
making it decay to zero in the long-time limit. From the 
m = mode results, displayed in panels (c) and (d), we 
see that the eigenvalues of the x- and the y-system are 
real and negative, excepting for the one closest to the 
origin of the Re(p) axis which is now exactly zero. Note 
that at T > 1/2, when m x n) = for all n e JV, the 
eigenvalue equations, ([53]) and (|54| . simplify for n > 2 to 



pm y n) 



Tn 2 m x n \ 
-Tn 2 m y n \ 



(55) 
(56) 
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while for n = 1, the equations are 

M m« = -Tm« + I Am (a)m«, 



/im 



^ = -Tm( 1 ) + iA m (a)mW. 



(57) 
(58) 



The system of equations, (f55]) - (|58|) is a set of independent 
equations, and may be solved quite easily to get the exact 
real eigenvalues [i. In particular, for X m (a) = 1 and T = 
1/2, one finds the eigenvalue /j, = 0. These eigenvalues 
agree with those shown in Fig. The zero eigenvalue 
makes the zero mode of fluctuations to neither grow nor 
decay at precisely the temperature T — 1/2. 

Based on results discussed in this section and in sec- 
tion QTJ we conclude that the linearly stable stationary 
state of the dynamics (0 in the large N limit is a state 
uniform in space but non-uniform in 6 when the temper- 
ature is less that 1/2. On tuning the temperature to a 
value larger than 1/2, the stationary state becomes uni- 
form in 6. Thus, our linear stability analysis in the limit 
TV —> oo allows us to predict that the system undergoes 
a transition at the temperature T = T C fl = 1/2. Indeed, 
simulation results shown in Fig. [5] for 7q* suggests a con- 
tinuous phase transition as a function of temperature, 
where we also show the theoretical curve for r?? obtained 



by combining Eq. 
the result 



30|) with Eq. ([33]) to give for T < 1/2 



h(m x /T) 



I (m x /T) 

for T > 1/2, on using Eq. (fTB")) . one has 



rS* = 0. 



(59) 



(60) 



dynamics in the continuum limit within a canonical en- 
semble of this so-called a-HMF model, by considering 
the Fokkcr-Planck equation for the evolution of the local 
density of particles. Our analysis shows that this model 
has a linearly stable stationary state which is always uni- 
form in space, being non-uniform in the internal degrees 
below a critical temperature and uniform above. Thus, 
the stationary state of the model is the same as that 
of its mean-field counterpart. We justified the mean-field 
dominance by performing linear stability analysis of both 
the uniform and non-uniform stationary solutions of the 
Fokker-Planck equation. 

The model we studied may be related to the Kuramoto 
model in the field of non-linear dynamical systems (IH - 
H^ |. Interpreting each of the N particles of our system 
to be a phase-only oscillator with phase 9i, the equation 
of motion ([7]) may be looked upon as the one governing 
the time evolution of the phase of the ith oscillator, in 
the presence of thermal noise. On considering a = in 
Eq. (0, the model in the absence of thermal noise is 
a particular limit of the Kuramoto model in which each 
oscillator has the same intrinsic frequency of oscillation. 
In the presence of thermal noise, Eq. ([7]) with a = 
is the generalization of the Kuramoto model considered 
in Ref. [l9j]. On considering a ^ 0, and attributing 
to each oscillator an intrinsic frequency sampled from a 
distribution, the model in the absence of thermal noise 
has been considered in [IH, [20|, Hl[ . In particular, a mean- 
field dominance similar to the one observed in this work 
has been seen in numerics in Ref. [15T ] . It remains an open 
problem to relate this dominance in these two scenarios. 



Exactly at T = 1/2, we have rg* = 0. 
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a=0.25,T=0.01 




0.0001 1 1 1 1 1 1 1 1 1 

5 10 15 20 25 30 35 40 

Time t 
a=0.5,T=0.05 




0.001 1 ' ' 1 1 1 1 1 1 

5 10 15 20 25 30 35 40 

Time t 



a=0.75,T=0.1 




0.001 1 ' ' 1 1 1 1 1 1 

5 10 15 20 25 30 35 40 

Time t 



FIG. 1: (Color online) Time evolution of the observables 
ro(t), ri(t), f*2 (i), and r?,{i) starting from an initial state 
{9i(0);i — 1,2, ...,7V}, where the 9i's are chosen uniformly 
in [— 7r,7r]. Thus, initially, the system is in the uniform 
state (|16[) . Here, the values of a and temperature T are 
(a) a = 0.25, T = 0.01, (b) a = 0.5, T = 0.05, and (c) 
a = 0.75, T = 0.1. From Table ((J), it follows then that 
the Fourier modes 0,1,2,3 in particular are linearly unsta- 
ble. Consequently, r*o (t) , r i (t ) , T2 (t ) , and (t) all grow in time 
from their initial values. However, the plots show that in the 
long-time limit, f*o(t) attains a value very close to unity, while 
Ti (t) , T2 (t) , T3 (t) all decay to a value very close to zero. The 
data in the plots are obtained from numerical simulations 
with N = 2 14 , and involve averaging over 100 independent 
initial conditions and dynamical realizations. The straight 
lines show the initial exponential growth with rates given by 
(T c , m — T), where the values of T^m's can be read off from 
Table [T] The agreement of the growth rates between theory 
and simulations is very good. 
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U°0=i.o 




3 , [^m(«)(1-T)-T] 1/2 ; , , } 

0.05 0.1 0.15 0.2 0.25 0.3 0.35 0.4 0.45 0.5 
Temperature T 



FIG. 2: (Color online) Plot showing as a function of T the 
values of m, that satisfy Eq. (|35[) , an d those that satisfy f|49[) 
at A m (a) = 1. The fact that the curves do not intersect at 
any T in the range < T < 1/2 shows that in this range, the 
value of m x that solves Eq. (|35|) at a given T does not satisfy 
Eq. 
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FIG. 3: (Color online) Real and imaginary parts of the eigen- 
values /j, of the x-system, Eq. (|53[) . and the y-system, Eq. 
(|54[) . for a m ^ mode such that A m (a) < 1 and for the 
m = mode such that A m (a) = 1. The temperature is 
T = 0.2. The parameter n max denotes the order of trun- 
cation of the eigenvalue equations. 
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FIG. 4: (Color online) Real and imaginary parts of the eigen- 
values [i of the x-system, Eq. (|53[) . and the y-system, Eq. 
(|54[) . for a m ^ mode such that A m (a) < 1 and for the 
m = mode such that A m (a) = 1. The temperature is 
T — 0.5. The parameter n ma x denotes the order of trun- 
cation of the eigenvalue equations. 
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FIG. 5: Plot showing rj* as a function of temperature. The 
data are obtained by starting from the uniform state (|16[) at a 
high temperature T = 1, and then changing adiabatically the 
temperature as a function of time t as T(t) = 1 — 7?£, where 
the rate R is taken to be such that Rr st <C 1, with r s t being 
the time to reach stationary state at a fixed temperature. The 
condition Rr s t <C 1 ensures that between two successive val- 
ues of T, the system gets enough time to attain stationarity. 
For a system of size N, taking r st = N 2 dt, where dt is the 
numerical integration timestep, we take R = 0.1/N 2 dt. 



